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Abstract. We study the conditions on the physical parameters in the Helfrich bending 
energy for hpid bilayer vesicles. The variation equation for embedded surface with a 
biconcave axisymmetric shape is analyzed in detail. This leads to simple conditions 
describing the solution and information about the geometry of the surface. 



Introduction 

In this article, a vesicle is represented by a closed surface S in M.^ with mean curvature 
H, surface area |S| and it encloses a volume V. Its geometric shape is modelled by 
minimizing a functional, sometimes called Helfrich functional, 



J^= [ {2H + co)^dS + X\J:\+pV 



with some physical constant parameters Cq, A, and p. The parameters carry the following 
meanings: Cq is the spontaneous curvature, A is the tensile stress, and p = — Pi is the 
osmotic pressure difference between the outer (po) and inner (pj) media. We take a sign 
convention that H is negative for the standard sphere. 

It has been observed experimentally long ago that a red blood cell is of biconcave- 
discoid shape (which will be defined mathematically later). And it is a quest to find the 
appropriate theoretical model for the bending energy. Historically, the early model of 
Canham, Q, is purely geometric and it is equivalent to the Willmore functional, |[W| , 



ch. 7], which equals with cq = A = p = 0, up to a constant. Certainly, from a differ- 
ential geometer's view, this cannot be the correct model because the unique minimum of 
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the Willmore functional for topologically spherical vesicles is the round sphere. This is 
also observed by physicists [pH|| . In fact, there are many important mathematical studies 



of the Willmore functional because of its geometric implications; for instance, the exis- 
tence of minimizers among a certain topological class by Simon and the conformal 



properties by Li and Yau [[LY|| . We expect that their works may contribute to a certain 
extent to a deeper theoretical understanding of the energy T . 

Helfrich takes physical condition together with the Gaussian curvature into account 



and proposes a modified bending energy, ||H1|| . The shape of blood cells and some other 



biological membranes is closely related to the formation of lipid bilayer vesicles in aque- 
ous medium (e.g. hquid crystal). The physical condition is based on the elasticity of lipid 
bilayer vesicles. According to the Gauss-Bonnet Theorem, the integral of the Gauss- 
ian curvature is a topological constant. Thus, within a certain topological class of E, 
Helfrich's bending energy can be reduced to T above. 

Many properties of T are yet to be discovered, though there are some experimental 



observations and numerical simulations, [|HDH] , |MB|| . The existence and uniqueness of 



its minimizer of a certain topology are still unknown. It is also not known whether the 
minimizer is symmetric in any sense. Answers to these questions require deep geometric 
analysis of the functional and studies in this direction are rare. Nevertheless, there are 
related works, such as, on similar functional, [S^, ^Y|, |NJ ; or on surface flows, |^ . 
The Euler-Lagrange equation corresponding to T is 



4Asif + 2(2if + co) 2{H^-K)-coH - 2XH + p = 
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In the past, much effort of physicists and biologists has been spent on studying axisym- 
metric solution to this variational equation, ||DH| , 0, |S^, PH2| , |0Y| , |N001| , |N002|] . By an 



axisymmetric surface, we mean an embedded surface S in which is rotationally sym- 
metric and has a reflection symmetry by the plane perpendicular to the rotational axis. 
With this additional assumption, the fourth order equation can be reduced to a second 
order one, usually referred to as the shape equation of axisymmetric vesicles. However, 
it is still unknown to scientists how the solution depends on the physical parameters 
and which parameters yield a solution corresponding to a biconcave surface. Our work 
is in this direction. We can derive conditions on the physical parameters for a solution 
having the biconcave shape. The conditions are easily expressed in terms of the cubic 
and quadratic polynomials denoted by 

Qit)=t' + 2cot^ + icl + X)t-^; 
R{t) = Q{t) - t\ 

Our main result can be stated as, 



Theorem. For any Cq, X, and p > such that every real root of Q is positive, there 
exists axisymmetric biconcave surfaces which satisfies the Euler- Lagrange equation of J-". 

It should be remarked that when cq > 0, A > 0, and p > 0, the condition is always 
satisfied. Helfrich's numercial simulation produces a biconcave shape resembling a blood 
cell when cq is positive. Furthermore, we are also able to numerically construct other 
interesting shapes when the condition is not satisfied. One is multiconcave and the other 
has no reflection symmetry. 
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This article is organized in the following way. In §1, we first give the differential 
equation for the revolving graph of an axisymmetric solution. The derivation of the 
equation is given in the appendix. Moreover, we formulate the problem of finding special 
solution corresponding to a biconcave shape surface. We also present several variants 
of the equation which will be useful later. The variants of the shape equation are then 
studied in §2 to show that our condition stated above is sufficient for getting the expected 
special solution for the problem. The analysis and the estimates of geometric quantities 
are discussed in detail here. Finally, since the solution and its reflection no longer form 
a graph at the reflection plane. We will show that the solution obtained in previous 
sections is still a solution across the reflection symmetry. It is sufficient to verify that 
it also satisfies the variation equation at the refiection plane. In the process, one more 
necessary geometric condition is obtained. 

Acknowledgement. We would like to thank our colleague K. S. Chou, who has been 
encouraging in our project and making valuable suggestions. 

1. The Equation, the Problem, and the Conditions 

An axisymmetric surface is a closed embedded surface E in with a rotation 
symmetry and a refiection symmetry by the plane perpendicular to the rotation axis. It 
is biconcave if there are exactly two components of negative Gaussian curvature. Without 
loss of generality, the rotational axis is labelled the 2;-axis and the plane of refiection is 
the xy-plane. Then the surface S can be obtained by revolving a radial curve about the 
2;-axis on the upper half plane and refiecting it to the lower half. Typically, a biconcave 
one is obtained by revolving and refiecting a curve shown in the picture below. 
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A cross-section of a biconcave axisymmetric surface 
(with Co = 1, A = 0.25, p=l) 

In other words, one may parametrize the upper part of E by X = (r cos 9, r sin 9, z{r)), 
with a function z{r) defined for r in some interval [0, r^], where z{r) > for r e [0, Too) 
and z{r^) = 0. There are natural boundary conditions imposed by the rotation and 
reflection symmetries of the surface. The obvious ones are z'{0) — and z'{r) — > — oo 
as r — > Too- However, more subtle ones arise from the regularity at the end-point Too, 
which we will discuss later in §3. For a biconcave surface, there is tm < such that 
z"{r) > if and only if r e [0, tm)- This is equivalent to require that z' has only a unique 
maximum and no other critical point. 

For a biconcave axisymmetric surface, the variation equation for the Helfrich func- 
tional is reduced to a second order ordinary differential equation, traditionally called the 
shape equation. With the notation iv{r) = z'{r), the shape equation is 

2r ,, 5rw ,2 2w' 

Ill — 111 — 

+ (l + w^/^ (l + ^/;2)V2 

210 + 2coW^ (cq + X)rw pr^ 

+ r(l + w2)3/2 + i + w^ ^ (1 + «;2)V2 ~ ~Y' 

We leave the derivation of this equation to the appendix to focus on the idea and analysis 
of the equation. 
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We are going to study special solution w{r) to the initial value problem on this 
variational equation with initial choice w{0) = and w'{0) = w'^ > 0. Specifically, we 
look for a solution with the additional requirements that 

CI: it is unimodal; and 

C2: there is a finite number r^r, such that w(r) - 



-oo as r ^ Too; and 



C3: -oo < / w(r)dr<0. 
Jo 



In [pH| , |JS|, |ZL|| , there are versions of the same equation written in term of the angle 
ip between the surface tangent and the plane perpendicular to rotational axis. However, 
it is convenient for our discussion to write the equation in the above form. 

With our notation of the polynomials Q and R, after multiplying with rw', the 
equation is equivalent to. 



where K(r) 



r 2/2 1 
r w 


/ 


_(l + w2)5/2_ 




r 2/2 1 
r w 


/ 


_(l + w2)5/2_ 





w 



(l + w;2)i/2 



-2 



+ r^w' R{K{r))] or 
+ rVQ(K(r)), 



w 



ryl + w'- 



-.. It will be seen that these groupings of the lower order terms 



are important in the analysis the equation. 

The understanding of n{r) also provides useful information about the solution. First, 
its derivatives are given by, 

K\r) 



w 



w 



r(l + w2)3/2 rWl + w^' 



K,"{r) 



w 



3ww 



l2 



2w' 



2w 



r(l + «;2)3/2 r(l + w;2)5/2 ^2(1 + ^2)3/2 



Then we have the equation for ^(r), 

—rn{rn' + k)2 



(2) 



2fl-r2/t2^ 



3/t' + 



rQ{K{r)) 



2(1 _ ^2^2^ 
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Let us end this section by remarking on a few geometric quantities in terms of w. 

w 

Firstly, n = — -^^=^= is the principle curvature in the meridinal (rotational) direction, 
rvl + w'^ 

w' 

Another principle curvature, the longitudinal one, is given by , , which occurs 

(1 + w^y^^ 

in K,'. Their product is the Gaussian curvature, which is expected to be positive at the 
rotational axis and reflection plane but negative somewhere along the circle defined by 
the zero of w. 



2. Analysis of the Equation 

In this section, we prove the conditions for the existence of the required special 
solution to the initial value problem on equation (0) described in the preceding section. 
The strategy is the following analysis on the equation. We study the principle curvature 
K, which is positive initially (at the rotational axis), i.e., k(0) = w'q > 0. If w'q is not too 
large in terms of the roots oi Q, k must decrease and eventually becomes zero at tq for 
some < ro < Too. After w becomes negative, it continues its descent and blows down 
to —oo at finite distance Too- This guarantees CI that w is unimodal and has a unique 
zero at vq. In order to verify that the solution satisfies the requirements C2 and C3, we 
establish, in terms of w'q, the estimates on tq, r^o and values of w, w' at these positions. 

Lemma 2.1. Let w'q > and R{t) < for t G [0,1^0]. If w > on an interval [0,ro), 
then K,{r) decreases on (0, tq). 

Proof. By continuity of k and that lim K,{r) = w'q > 0, there exists e > such that 
i?(/t(r)) < 0, w > 0, and w' > on (0, e). Integration of the equation (^a) on (0, e) gives 

2/2 2 

r w w 

-, < , on (0,£). 

(l + ti;2)5/2 VTT^ 
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w' W 



Thus, n'ir) = — , < for r G (0,e). This also implies that /t(r) 

remains in the interval (0, w'q] for r G (0, e) and hence the argument works before we hit 
the first zero tm of w', i.e., the first critical point of w. This shows that ri = inf{r G 
(0,ro) : K'(r) > 0} > Vm > if the set is nonempty. 

Now by the smoothness of n, we conclude that ^(ri) G [0,Wq], n'{ri) = and 
> 0. Putting this into the equation (|^), we have 

° - "^'^ ^"^^ - 2(l-r,V(,,)) < 0- 
This is a contradiction and hence the set {r G (0,ro) : ^'(r) > 0} must be empty. □ 

By the same token, we can also establish a criterion for a growing solution. Since 
we don't need it for our further discussion, we will omit the proof. 

Proposition 2.2. If cq > 0, p > 0, Wq > and R{w'q) > 0, then w is increasing and 
blows-up to +00. 

For our future discussion, we denote a few quantities which depends only on the 
polynomial Q and the initial data w'q as follow: 

I^i^Wq) = - min {Q{t) : < t < Wq} = max {-Q{t) : < t < Wq} ; 
5+{wq) = mm{-Q{t) : < t < w'q} = - m&x{Q{t) : < t < Wq} ; 
6_ = min {-Q{t) : t < 0} = -max{Q{t) : t < 0} . 

Note that lim S+{w'q) = p/2 = lim /u(wo), and 6- is independent of w'q. 

We will now show that under reasonable condition, the solution will not blow up to 
+00. It is because the cubic lower terms are dominated by S+. 
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Lemma 2.3. Suppose that all real roots of Q are positive. If Q < on [0,Wq], i.e., w'q 
less than the smallest real root of Q, then there is ro > with Tq < such that on the 
interval (0, Tq), 

w > 0, w{ro) = 0, and w'(ro) < — -^^^o- 

8 

In addition, if 6Awq^ < 275+, then 

Proof. Consider the principal curvature ^(r) = — j^=^= and rewrite equation (|l]) or 
(0a) as 



= n ^ 2^5/2 + '^ O • 



By lemma |2.1| , ^'(r) < on a neighborhood of and hence K{r) < w'q near 0. Using the 



assumption on Q and the positivity of w near 0, we have 



Therefore, for small r, 



(3) K'(r) < -^r^, 



8 



and 



(4) n{r) <w',- ^r\ 

These inequalities show that k remains in (0,Wq) as long as w > 0. In turns, they 

themselves hold as long as w > 0. Clearly, they guarantee the existence of ro with 

Tg < ° such that K(ro) = and so w{rQ) = 0. Then by (|), we have w'iro) < — 5+rQ/8. 
(>+ 

This proves the first statement of the theorem. 



10 
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To prove the second statement, we temporarily let a = w'q and h = 5+/16. Then on 
[0,ro), 



W 2 

< a — or . 



r\/l + 

With the assumption that 4a^ < 27b, it can be written as 

r(a — br"^) 

w < 



1 — r^fa — br'^)'^ 



Therefore, 



ria — br"^] 

w < 



276 



and 



war < I — , dr 



JO .11 _ M 



< 



1 Aw'^ QAw'J'V^'^ 



27 



This completes the proof of the lemma. □ 
Immediately from the lemma we have the following 

Corollary 2.4. Suppose that all real roots of Q are positive. Then 

K.2 



lim sup — 7 < — and lim sup — ^ [ ^'^ dr 

w'^-^o Wq p ^^^0 w'q Jo 



< -. 

p 



In the above, we obtained an upper bound for ro, where w first hits zero. Next, a 
lower bound is established, which is essential for future estimates. We note that from 
lemma |2.3|, w has at least one maximum in the interval (0,ro). Since w is increasing at 



the beginning, the first critical point must be a maximum point. From now on, we let 
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Lemma 2.5. Let tm G (0, tq) he the first critical point of w. Then 



liminf ^ > liminf ^ > ^ 



!m >^ 



liminf^>-2. 



Proof. Using equation (^), k > on [0,ro), and Tq < 16w'q/6+, one can check that 



1 — r K{r) >^ on(0,ro). 



Putting together the assumption on Q, the fact that k < Wq on (0,ro), and the inequal- 
ity (H) into the equations (@), one obtains 

{rK)" > —^rK,{rK + nf + - ^r. 
z4 o Zq^ 



2i 



Note that in the above, /i = max {— (^(t) : < t < Wg} > 5+. 



Let y = (r/t)' 



:i + w2)3/2 



and A' = i^^^ > iLJl^ > ^ > 0. We then 



have 



-Wn, 

> ^r. 



Wq 



Aw'o 



Integrating from 0, it yields 



(5) 



y 



W'n \ Aw'n 



arctan — > arctan — — 



.AJ - \A J 4^ 

Taking r = rj\,/ and multiplying by A/w'q, as y{rM) = and tq > rjv/, we have 



A 



A^l > A^rl > 4e ■ 4 arctan [^)>4^1 



/2 ■ 



W, 



A 



3^2 



Then, 
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Considering the limiting situation, by lim fi = lim 5+ = p/2 and lim ^ = 1, the first 



result follows. 



For the second estimate, we first take r ^ ro in apply Tq < 16wq/(5+ to obtain 



Ajw'jro) - w'o) 
+ w'{rQ)wQ 



> tan 



Then, 



w'q 



cos 



> 



> 



-4Aw: 



i2 



hi 



AAw. 



COS 



hi 



/3 



This gives the second estimate immediately by letting — ^ 0. 



Combining the previous lemmas, we have the following 



Corollary 2.6. For any e > Q, there exists (3 > Q such that if w'q G {0,j3), then 



32 



3p Wq Wq p 



32 



and 



w'(ro) 2 

Wn 3 



□ 



The rM that we discussed in the previous lemma is in fact the unique critical point 
of w in (0, ro) provided w'^ is sufficiently small. 



Lemma 2.7. Under conditions of lemmas p.3| and p.5| , if w'q is sufficiently small, then 
w has a unique maximum at rM G [0, rg) and has no other critical point. 



Proof. Suppose w has another critical points in (0,ro), then according to the definition 
of rM, w attains a positive nonmaximum critical at G {rM,ro), i.e., w'{rm) = 
and w"{rm) > 0. By the corollary |2l6| , given any £ > 0, > rj^j > (|f ~ ^) ""^o ^'^^ 
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sufficiently small w'q. On the other hand, substituting w'{rm) = and w"{rm) > into 
the equation ([^), we have 

2Wm 



< 



+ Q{K{r,n)) 



where Wm = u){rm)- Since < ^(rm) < w'q, thus, for sufficiently small w'q, we also have 
Q{K{rm)) < — p/2 + e. This leads to 

2 p 
< ^ ^ + £ 



f 2 



which is clearly a contradiction since e is arbitrary. 



□ 



Next, we claim that under our condition on Q, after hitting 0, w decreases and goes 
to negative infinity in finite distance. 



Theorem 2.8. Suppose that all real roots of Q{t) are positive and ro > is given by 
lemma p.3| . Then w' < for r > tq and there is a finite number Too > tq such that 
— oo. Moreover, we have 



lim w{r 



roo-rQ< 



IT 



2^Sri\w'iro)\ 

where 6 = min and (5_ = min {—Q{t) : t < 0}. 



Proof. For convenience, let us temporarily denote v = —w for r > rg. Since K,{r) < and 
f ' > on (ro,ro + e) for some £ > 0, we have —v'Q{K,{r)) > 6^v' on the same interval 
and the equation (|I|b) gives 

' -2 



(6) 





/ 


-2 


_(l + t;2)5/2^ 




_(l + i;2)V2_ 



- rVQ(/t(r)) > 



(1 + ^2)1/2 



+ S-r-^v 



3„,/ 
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on (ro, ro + e). In particular, 



2/2 

r V 



(1+^2)5/2 



> 



(l+i;2)V2 



which gives 



roV (ro) > 



(1+^2)5/2 V u; _ (1^^2)1/2 

This imphes that v' does not vanish and therefore 



+ 2 > 0. 



(7) 



rv' > 



rv 



> rov'{ro) > 0. 



(1 + t;2)5/4 

This, in turns, imphes that v is increasing, k remains negative and @ holds as long as 
V is defined. Substitute into the equation (Q) again, it gives 



2/2 "I ' 

r V 



> 6-r v' > (5-rof'(ro)r > d-r^v' {rQ)r 



(1 +t;2)5/2_ 

for all r > ro such that v is defined. Now, integrating the above from ro, we conclude 
that 



2 /2 

r V 



>roV(ro)^ + -^-^°"'^^°^-^ 



r - rn 



(1 + ^2)5/2 - V - ■ 2 

In the case that 5+ > 45_, we have f'(ro) > (5_rQ/2 by lemma p73 

5 



(l + t;2)5/2 



> row'(ro) t;'(ro) - — rg + 



> 5-^jt/(ro)^2 



5_rgf'(ro) 



Otherwise, for ^ = ^ < 1, one has rQt>'(ro)^ > \C^~v'{rQ)rQ. It follows that 



^2^/2 



(1+^2)5/2 



\ 2 1/ \2 I ^-'^qV (ro) ,2 2\ 

>rof (ro) + ^ (r - rg) 

2 // ^2 , C^~rlv'{rQ) 2 2N , (1 - C)'^-'"0^'('^0) . 2 2N 

= rot;(ro) + (r - rg) + (r - rg) 



> r^v (ro) 



i/(ro) 



C5-t;'(ro) 
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Therefore, 



(l + t;2)2 - (1 +^2)5/2 - 2 

Hence, with 5 = min |<^-} , one has 

/ 

> \j5rlv\rQ) 



v' ^ v' 



l + v^ - (l + t;2)5/4 
as long as v is defined. Further integrating from rg gives 



arctanf > \j6r^v'{ro){r — tq) 



This clearly shows that v blows up before r goes to infinity. Furthermore, the upper 
bound of Too follows. □ 



Corollary 2.9. Under the same conditions of theorem |2.8|, we have, 



oo> I \w{r)\dr > B ^\ogl —-^ rr ) > ^(^oo - ro)^ 

Vcos(5(roo - ro))/ 2 



ro 



where B = ^6r^v'{ro) > ^^r^ > -l^Srl 



Proof. From (^, by simply dropping the term —r^v'Q{K,{r)) > 0, we have 

r\'^ ( 1 \ 



> 2 1 



(1 + 1;2)5/2 - \^ yiT^ 
After multiplying by vT+^(vT+^+ 1), it becomes 



(Vl +t;2 + 1) / ^ > 2(Vl + t;2 + l)(Vl + t;2- 1) = 2t;2. 
(1 + f 

Moreover, 2\/l + t>2 > (vTT^+ 1), so one has 



/2 

1 + ^)3/2 



2 V ^2 



Taking square root and integrating, it yields 



v' 

dr > I vdr 



16 
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in which the left hand side is obviously convergent. The lower bound follows easily from 
the theorem. □ 

It can be seen that the area of negative part has a lower bound of order Wo(^oo — ''"o)^- 
We wish to estabhsh that this dominates the area of the positive part. Let us consider 
the equation of v{r) = ~w{r) for r > tq again, i.e., 



2/2 

r V 



(1 + ^2)5/2 



< 



< 



(1 + ^2)1/2 

2 1 ' 

(l + t;2)i/2 

(l + t;2)l/2_ 



+ 



2cof ^fV (cg + X)r'^vv' pr^v' 
1 + ^ (l + t;2)i/2 

1 + (l + t;2)i/2 + 2 



+ 



2\co\r^ + {cl + \)rl + P'"^ 



I 

V . 



Therefore, after integrating from tq. 



2 /2 

r V 



< TqV (ro) + 



(l + ^;2)5/2 - V (1+^2)1/2 



+ 



2\c,\r^ + {cl + \)rl + ^-^ 



l + 2\co\r^ + {cl + X)rl + ^ 



2„,/2 



(l + t;2)5/2 



l + 2|co|roo + (c^ + A)r^ + ^r^ 



From this, we have proved a comparison between the growth of w' and \w\. 



Proposition 2.10. Ifw{r) is a solution to equation (|T]) which blows down to — oo at Too, 

/2 

7 7 7 7 

then - — ^rTTTTT bounded. 



\w\ (1 +M;2)5/2 



Continue with the analysis, after integration and applying Holder inequality, we have 



2 ^ 2 



1 <r, 



ro (l+t;2)5/4 y - % (l+t;2)5/2i,„ 



,/2 



< rlv'{ro)\r^ - ro)' + (r^, - Tq) 1 + 2|co|roo + (eg + A)r^ + ^r^ 
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Letting x = r^o — and using ^]6| that Tq < f^Wg ^^'^ "^^'(^0) < 2wo, we have 



2 ^ 64 ,3 2, 



2 P S 

Oo + aix + + -x 



'"0 



where 



ao = l + 2|co|ro + (c2 + A)r2 + |r^, 

ai = 2|co| + 2(co + A)ro + — rg and 

/ 2 XX 3p 
^2 = (cq + A) + — ro- 



We now claim that there exists C > independent of Wq such that f^^ v > Cw'q for 
sufficiently small w'q. Otherwise, for any e > 0, there is a < e such that 



SWq > I V. 



Then by corollary |2]^, we have 



ewr. > — i=rnX . 



Using the lower estimate of Tq in terms of w'q, we conclude that 



2 2V2w' 2V26+ 
X < e— ^ < —^e. 



Putting this into (H), we have 



ciu^o < C2w'o'e + C3£^/' {1 + o{w'o^^^) + [2|co| + o{w'q'^^) 



+ 



.2. , E^3^3/2 



for some positive constants Ci, C2 and C3. This is impossible since this implies 



< ci < C2w'Q^e + c^e 



_3/2 



1 + o(e^/^) + o(w, 



which is a contradiction. 
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Theorem 2.11. With the conditions of the theorems |2.3| , |2.5| , and |2.8| , we have 



-oo < / w(r) dr < 
Jo 



for w'q sufficiently small. 



3. A Necessary Geometric Condition 

The fact that the surface has a reflection symmetry through the xy-p\a.ne actually 
provides some geometric conditions on the surface. Suppose we have a function z{r) with 
w = z' = Zr satisfying the equation (|l|), we must check that the surface obtained by 
revolving the function —z{r) is also stationary to the functional JF. For this purpose, we 
may parametrize the surface by an even function r{z) with z G [—a, a], namely, 

'X(5,z) = {r{z) cose, r{z) sine, z), {e,z) G [0, 27r] x [-a, a], 

r{z)=r{-z), r{0)=r^. 

In this section, we will use Zr or Wr to stand for derivatives wrt r, i.e., z' and w' previously 
in order to distinguish from derivatives wrt z denoted by another subscript. 

The Helfrich functional for r(z) and —r(z) near r^o can be written as, 



— = lim ,^ 

27r e^Oj^a 



La + Is ( [^^^ + """^'l + f ) 



We then study the critical function of the functional. Assume that the variation 
6{r{z)) = ip{z), then the above is written as 



27r 



lim j {F2(pzz + Fiifz + Foif) dz; 
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where (with similar calculations as those in the appendix) 
2rr^^ 2 , 2cor 



2 



(r2 + 1)5/2 (^2 + 1)3/2 ^2^1' 



^ ~ (r2 + 1)7/2 ^(^2 + 1)3/2 + (^2 + 1)1/2 + (^2 + 1)5/2 (^2 + 1)2 ' 
„2 



Since F2 and Fq are even and Fi is odd, after integrating by parts, we have 

^ = r [(F2),,-(Fi). + Fo]v.dz-lim [-(F2),(e) + Fi(5)] [^,(5) - ^,(-e)] . 

The function r{z) here is the inverse function of the function z{r) in the previous sections. 
Thus, Tz and r^^ above can be written in terms of w by chain rule. Since w{r) = z'{r) is 
a solution to the equation (|1|), we obtain 

{F2),, - (Fi), + Fo = 0. 

def 

Furthermore, we may denote //(r) —(^2)2 + Fi. Then, in terms of w, 



99 9 
V = M , .■ 9N.^/9 1. . o - 2cow - (cq + A)rVl + + ^ 



(l+y;2)5/2 ^yi" 

Consequently, 5(J^) = if and only if //(r) is bounded in a neighborhood of Too. Note 
that 



r]r -Wl W 2CoWr , 2 , ,^ Vl + 

= ~ = + 0/1, o 7^ 2Y ~ ^0 + A +pr 

w ^(l + w'')^/^ r^vl+w^ w(l + w^) w 

is bounded according to the analysis in the previous section, specifically by proposi- 
tion |2.1(]| . Therefore, ?7(r) is also bounded and hence, we have established the following 



theorem. 

Theorem 3.1. Let z = z{r) > be a function on [0,roo] with z{r) = if and only if 
r = Too- If w = z' is a solution to the equation then the surface obtained by revolving 
the curves z{r) and —z{r) is stationary to the functional T . 
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Furthermore, according to this boundary condition, by considering hm 
0, there is a requirement on the geometry of the surface. 



r]{r) 



Proposition 3.2. Let E be an axisymmetric stationary surface of T. Ifr^ is the radius 
of the circle of intersection of E and the reflection plane, then the Gaussian curvature of 
every point along this circle is given by K{r ooY = — Q ( — ^ • 



4. Appendix: Derivation of the Main Equation 

Let the surface of revolution E be parametrized by X(r, ^) = (r cos ^, r sin ^, 2;(r)) 
for (r, ^) e [0, Too] X [0, 27r]. Then, with w{r) — z'{r), the mean curvature H is given by 



H 



/I la' 
w H — w H — w 



2{l + w^) 

The area element is 27rr\^T+^ dr. With the assumption that z{roo) — 0, after inte- 



grating by parts, the volume enclosed by E is — tt / r w(r) dr. Therefore, to find the 

^0 

critical point of we may consider the variation on 



27r JO 



(2H + cq? + X] rVT+w^ dr - ^ 1'^ r'^wdr. 
J 2 Jo 



Let 5 be the variational operator and 6{w) = ip where ip{r) is a smooth function with 
compact support. Then ST/2'k — I1+I2+I3 where 



Ti= / 2(2H + Co) 5(H) rVl + dr; 
Jo 

I2 = \{2H + Co)' + a] rS{Vl+^) 

Jo L ^ 

' \2H + Co)' + A 



dr 



T3 = -y^ r5Hdr^-j^ ^ 



rwcp , 
r , dr; 

r''(f dr. 
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We also have 



6{2H) 



(1+^2)5/2 
1 

(1+^2)5/2 



(1 + w'^)ip' H — (1 + w'^)ip — 3ww'(p 
r 

(1 + w'^)w" - 2,ww'^ + -(1 + w'^)w' - ^w{l + w'^f 



Therefore, 



2r{2H + Co) 
(1 + ^2)2 

2r{2H 



(1 + w'^)ip' H — (1 + w'^)ip — 3ww'(p 
r 



Jo dr 



Co 



Jo 



"•oo 6rww'(2H + Co) 



(l + «;2) ^ 



(1 + W"')^ 

-2r{2Hy 2rww'{2H + cq) 



—2rw" 



(1 + w'^y 



+ 



Qrww 



,2 



ipdr 
2w' 



+ 



2w 



(l + w2)5/2 (1+^2)7/2 (1+^2)5/2 ^(1 + ^2)3/2 

2corww' 2rww''^ 2w'^w' 



(1 + ^2)2 (1 + ^2)7/2 (1 + ^2)5/2 



ipdr. 



The quantity I2 is expanded into the following. 



I2 



J 

Jo 

f 

Jo 



[{2H + co)^ + X\rw 



(1 + «;2)V2 
rw 



ip dr 



(1 + w2)7/2 

2cor'u; 



+ 



(1 + W2)2 



rww 



/2 



+ 



2w'^w' 



+ 



w 



(1+^2)7/2 (l + u,2)5/2 r(l+w2)3/2 

2cQrww' 2coU'2 (cq + X)rw 



+ 



(1+^2)2 (l + w;2) ^/r^p^ 



ipdr. 



ip dr 
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Summing the above quantities, we have 



ZTT 



rr 

Jo 



-2rw" 



hrwW 



,/2 



2w' 



+ 



2w + w'^ 



(l+«;2)5/2 (1 + m;2)7/2 (1 + ^2)5/2 r(l + ^/;2)3/2 
,2 /'„2 I ^^2" 



^ 2cow^ ^ (cq + A)rw pr 



if dr. 



(1 + W2)V2 2 

Therefore, the variational equation is given by (integrand) + constant = 0. With the 
assumption that w{Q) — 0, taking r — > 0, one gets 

2w 

—2w'{Q) + hm h constant = 0. 



This concludes the derivation of the variational equation. 
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